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ABSTRACT 

Background: Community-acquired methicillin-resistant Staphylococcus aureus 
(CA-MRSA), a novel strain of MRSA, has recently emerged and rapidly spread 
in the community. Invasion into the hospital setting with replacement of 
the hospital-acquired MRSA (HA-MRSA) has also been documented. Co- 
colonization with both CA-MRSA and HA-MRSA would have important clinical 
implications given differences in antimicrobial susceptibility profiles and the po- 
tential for exchange of genetic information. 

Methods: A deterministic mathematical model was developed to characterize 
the transmission dynamics of HA-MRSA and CA-MRSA in the hospital setting 
and to quantify the emergence of co-colonization with both strains. 
Results: The model analysis shows that the state of co-colonization becomes 
endemic over time and that there is no competitive exclusion of either strain. 
Increasing the length of stay or rate of hospital entry among patients colonized 
with CA-MRSA leads to a rapid increase in the co-colonized state. Compared 
to MRSA decolonization strategy, improving hand hygiene compliance has the 
greatest impact on decreasing the prevalence of HA-MRSA, CA-MRSA and the 
co-colonized state. 

Conclusions: The model predicts that with the expanding community reservoir 
of CA-MRSA, the majority of hospitalized patients will become colonized with 
both CA-MRSA and HA-MRSA. 
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Introduction 

Methicillin-resistant Staphylococcus aureus (MRSA) has been traditionally considered 
a hospital-acquired bacteria and is implicated in the great majority of infections acquired 
in the hospital [ij]. The documentation of a novel MRSA strain, which has emerged in the 



community and has subsequently spread into the hospital, 



las led to a re-evaluation of the 



transmission dynamics of MRSA in the healthcare setting . Several population-based 
surveillance studies have documented that the community- acquired MRSA (CA-MRSA) 
may be replacing the hospital-acquired MRSA (HA-MRSA) 0,151,0]. Mathematical 
models corroborate these findings and predict that there will be competitive exclusion 
of HA-MRSA strains by CA-MRSA over time 7|. Previous studies have assumed that 
individuals can only be colonized or infected with either HA-MRSA or CA-MRSA. Data 
suggest however, that individual colonization with multiple Staphylococcus aureus strains 
occurs 8(1 . CojColonization with multiple strains of other bacterial species has also been 
documented [9]. Understanding the emergence and spread of co-colonization with both 
CA-MRSA and HA-MRSA would have important clinical implications given differences in 



antimicrobial susceptibility profiles and virulence properties between these two strains [10|. 
Co-colonization can also result in the horizontal transfer of mobile genetic elements between 
strains, such as antimicrobial-resistance or virulence determinants. These events may 
lead to the emergence of MRSA strains with novel biological properties. A mathematical 
model was developed to understand the emergence and spread of co-colonization with 
both CA-MRSA and HA-MRSA among hospitalized individuals. This model extends a 
previous model characterizing the transmission dynamics of CA-MRSA into the hospital 
setting, which assumed that patients could only be colonized or infected with either CA- 
or HA-MRSA [7|. Key factors, which contribute to the spread of antibiotic-resistant 
bacteria and their impact on the emergence of a co-colonized state in the hospital setting, 
were analyzed through model simulations. The impact of an increased influx of patients 



4 



harboring CA-MRSA into the hospital was quantified using data from population-based 
surveillance studies, which document an expanding community reservoir of CA-MRSA. 
Differences in length of stay (LOS) among patients harboring CA-MRSA were also analyzed 
since patients infected with CA-MRSA can present, with severe infections leading to longer 
LOS. Lastly, infection control strategies aimed at hmiting the spread of MRS A and their 
effect on the emergence of the co-colonized state were evaluated. 

Methods 

Mathematical Model 

A deterministic model was developed to characterize the transmission dynamics 
of HA-MRSA and CA-MRSA in the hospital setting and to quantify the emergence of 
co-colonization with both CA-MRSA and HA-MRSA among hospitalized patients. 

Individuals in the hospital are in four exclusive states: susceptible (S), colonized with 
either CA-MRSA (C), HA-MRSA (H) or both CA- and HA-MRSA (B). The infected state 
is not included. Patients enter the hospital as -S, C or H. To understand the emergence 
of the co-colonized state, the model assumes that there is no co-colonization at baseline 
and that patients do not enter the hospital in the B state. Patients leave the hospital 
via death or hospital discharge in all four states. Through contact with a contaminated 
healthcare workers (HCW), susceptible patients becomes colonized with either CA-MRSA 
at a rate (1 — r])l3c or HA-MRSA at a rate of rate (1 — r])(3H- Here, rj signifies compliance 
with hand hygiene measures with 0<77<1,77 = corresponds to no compliance and 
rj — 1 corresponds to perfect comphance. Transmission rates of CA-MRSA and HA-MRSA 
are given by f3c and (3h respectively. Once in the C ot H state, patients can transition 
to the co-colonization state, B, through contact with a HCW, contaminated with either 
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HA-MRSA or CA-MRSA, with rates (1 - r])pcH or (1 - r])pHC, respectively (see figure [T] 
and appendix for model details). 

The LOS among CA-MRSA colonized patients was assumed to equal to the LOS of 
susceptible patients in the baseline model. The LOS for the co-colonized compartment 
starts after acquisition of the second strain and, as a simplification, is set equal to 
the larger of the LOS for patients colonized with CA-MRSA or the LOS for patients 
colonized with HA-MRSA. Since patients colonized with MRSA can develop an infection 
during their hospitalization, thereby substantially prolonging their LOS, simulations were 
performed to quantify the impact of an increasing LOS on the transmission dynamics of 
HA-MRSA, CA-MRSA and the emergence of the co-colonization state. Simulations were 
also performed to determine the impact of an increase in the percent of patients entering 
the hospital already colonized with CA-MRSA. Model simulations evaluating the impact 
of hand-hygiene and decolonization of MRSA colonized patients, two key infection control 
strategies, were also performed. Since several different decolonization strategies are available 
with varying efficacy, the decolonization parameters of patients colonized with CA-MRSA 
(ac), HA-MRSA {an), or co-colonized (a^) ranged from 0% efficacy (no decolonization) to 
100% efficacy. 

The mathematical model assumes that the likelihood of transmission of CA-MRSA and 
HA-MRSA are equal [jSc = Ph)- In vitro data suggest that the growth rate of CA-MRSA is 
faster than that of HA-MRSA for certain CA-MRSA strains. This biological difference may 
allow CA-MRSA to have an advantage towards colonization and subsequent transmission 
compared to HA-MRSA. To understand the implications of a greater transmission potential 
among CA-MRSA strains, the baseline model and above simulations were re-analyzed with 

(3c>(3h flU. 

Parameter estimates were obtained from infection control data, microbiology data, and 
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patient data from a 400-bed tertiary care hospital with approximately 25,000 admissions 



per year. Estimates were also obtained from pub 
of CA-MRSA or HA-MRSA (Table 1) Q, Q, 12, 



ished studies focusing on the epidemiology 

3- 
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15, 



Results 

Transmission dynamics of HA-MRSA and CA-MRSA and co-colonization 
Baseline model 

A baseline model was developed to quantify the prevalence of patients colonized 
with HA-MRSA, CA-MRSA, or both strains over time. To understand the underlying 
transmission dynamics of HA- and CA-MRSA and the emergence of co-colonization with 
both strains, this baseline model assumes that there is no entry of patients who are already 
colonized with MRSA into the hospital (all entering patients are susceptible). Analysis 
shows that when the basic reproduction numbers satisfy RO^ > RQ^ > 1 (see appendix), 
the state of co-colonization with both HA-MRSA and CA-MRSA becomes endemic over 
time. Assuming P = Pc = Ph = Phc = Pch and a = ac = c^h = c^b, the model analysis 
also demonstrates that there is no competitive exclusion of either strain, when both RO'" > 1 
and -RO^ > 1. Over time, the prevalence of patients colonized with HA-MRSA exceeds the 
prevalence of patients colonized with CA-MRSA since RO^ > RO^. The greater RO^ value 
for HA-MRSA compared to the RO'-' value for CA-MRSA reflects the longer LOS among 
HA-MRSA patients, which results in greater opportunities for HA-MRSA transmission 
(figure E]). 
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Simulation 1: increased transmission of CA-MRSA and HA-MRSA 

Increasing patient-to-patient transmission through patient contact with HCW 
contaminated with either HA-MRSA or CA-MRSA results in a substantial increase in the 
percent of co-colonized patients. As transmission increases, more patients are colonized 
with MRSA and less are susceptible, and therefore more patients that are colonized with 
individual strains become co-colonized. Above a threshold value of /3, the percent of 
patients co-colonized with both strains exceeds those colonized with either CA-MRSA or 
HA-MRSA (figure [3D. 

Simulation 2: increasing the influx of patients colonized with CA-MRSA into 

the hospital and their LOS 

Increasing the rate of admission of patients colonized with CA-MRSA or increasing 
their LOS results in an increase in the co-colonized state (figures H] and [5]) . LOS has 
a substantially greater impact on the prevalence of co-colonization compared to rate of 
admission. Even a minimal increase in LOS from the baseline value of 5 days to 8 days 
leads to the majority of colonized patients represented by the co-colonized state. The 
predominance of the co-colonization state when LOS increases is explained by the increased 
opportunity for HA-MRSA colonized patients to become colonized with CA-MRSA and 
become colonized with both strains. 

Simulation 3: interventions to prevent transmission 

The effect of two standard interventions aimed at preventing the transmission of MRSA, 
improving compliance with hand-hygiene and maximizing the efficacy of decolonization of 
patients with MRSA, were evaluated in simulations which included an inffux of patients 
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harboring HA-MRSA or CA-MRSA. Both interventions decrease the percentage of colonized 
patients in all three states. Compared to decolonization, improving hand-hygiene has the 
greater impact with even small increases in compliance having a substantial effect in the 
overall prevalence of MRSA. Since there is a constant influx of colonized patients into the 
hospital setting, CA- and HA-MRSA are never eradicated from the hospital even at 100% 
hand-hygiene compliance or 100% decolonization efficacy. In contrast, the absence of an 
influx of patients who are co-colonized into the hospital explains the extinction of the B 
state when these two interventions are at 100% compliance or efficacy. As hand-hygiene 
compliance increases, the total percentage of patients colonized decreases monotonically. 
However, simulations show that after 2 years, the percentage of patients who only have 
HA-MRSA increases until hand-hygiene reaches about 45% (figure [6]). 

The explanation for this finding is as follows. Susceptible patients move to the single 
colonized state with a rate {1 — r])j3 and from the single colonized state to the both state 
with a rate (1 — r^)/?. Therefore there is a quadratic effect of hygiene on transmission 
to the co-colonized state. When hand-hygiene compliance is low (and the co-colonized 
compartment is large), the quadratic effect of increasing hand-hygiene compliance, causes a 
rapid reduction of the co-colonized compartment, and increases the population of susceptible 
patients in the equilibrium. Transmission to the single colonization compartments {C and 
H) is dependent both on r] and S by the term ^^-^^S. Until hand-hygiene reaches about 
45%, the rise in 5* due to the quadratic effect on the co-colonization state is stronger than 
the reduction in transmission due to increasing rj. Therefore, the sizes of the C and H 
compartments increase. 
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Simulations assuming greater transmission potential for CA-MRSA compared 

to HA-MRSA 

The overall results of simulations with Pc = 1-33/5// were similar to those with Pc = Ph 
except for a greater and more rapid increase in CA-MRSA and co-colonized patients (data 
not shown). 

Discussion 

The transmission dynamics of MRSA in the hospital setting are complex and require 
the analysis of numerous interrelated and dynamic factors. The emergence of CA-MRSA 
and its invasion into the hospital setting has led to further complexities in understanding 
not only the spread of MRSA, but also the selection of effective antimicrobial therapies and 
preventive strategies. Since epidemiological studies cannot fully address these complexities, 
a mathematical model was developed to specifically quantify the emergence of individuals 
co-colonized with both CA-MRSA and HA-MRSA strains. 

The deterministic model shows that over time, there will be no competitive exclusion 
of either strain but that both CA-MRSA and HA-MRSA will co-exist in the hospital 
setting. The model also shows that individuals co-colonized with both strains will increase 
in prevalence over time and will predominate over individuals who are colonized with 
either CA-MRSA or HA-MRSA. These findings have important implications. First, clinical 
cultures will usually identify only one strain, either CA-MRSA or HA-MRSA. Mathematical 
models have shown that 20 colonies_per patient need to be sampled to reliably estimate 
the occurrence of multiple strains [17,]. Since sampling of multiple colonies is not feasible, 
co-colonization and polyclonal infections will therefore not be routinely detected. The 
different antimicrobial susceptibility profiles of CA-MRSA and HA-MRSA and identification 
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of only one strain may therefore lead to incorrect antimicrobial therapy. The emergence of 
multidrug-resistant strains of CA-MRSA, with susceptibility profiles similar to HA-MRSA 
may however minimize difficulties in selection of the appropriate antimicrobial among 
co-colonized patients [l8|. 



Factors which increase the reservoir of CA-MRSA in the hospital setting were shown 
to have substantial effects on the magnitude of the co-colonized patient population. 
Increasing the influx of patients harboring CA-MRSA into the hospital and increasing 
their LOS resulted in a rapid increase in the number of co-colonized patients. Previous 
mathematical models have also shown that these two factors are central to the dissemination 
of antimicrobial-resistant bacteria P, [l9]. Our model shows that even small increases in LOS 
of only few days from a conservative baseline estimate of 5 days among CA-MRSA patients 
resulted in the rapid predominance of patients colonized with both strains. CA-MRSA 
is associated with severe infections and several studies have documented that CA-MRSA 
has become the predominant MRSA strain implicated in the great majority of nosocomial 
blood stream infections and surgical site infections 4 , Q, [2lj] . These nosocomial infections 
would therefore substantially extend the LOS of patients harboring CA-MRSA. 

Our model revealed the effects of two interventions targeting the prevention of MRSA 
spread: improving compliance with hand-hygiene and decolonization strategies. Both 
interventions decreased the prevalence of HA-MRSA, CA-MRSA and co-colonization. As 
shown in previous models, improving hand-hygiene compliance had the greatest effect with 
only small improvements in compliance Q]- 

Several assumptions were made to simplify the model. First, antibiotic pressure and 
its effect on the transmission dynamics of CA-MRSA and emergence of co-colonization 
were not assessed. Given the different susceptibility profiles between HA-MRSA and 
CA-MRSA, selective antibiotic pressure may alter the transmission dynamics between 
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HA-MRSA and CA-MRSA. Second, environmental contamination was not included since 
there is a paucity of data regarding differences in contamination of inanimate surfaces 
between CA-MRSA and HA-MRSA. Future models will need to include these factors. 
Our main model assumed that the likelihood of transmission between CA-MRSA and 
HA-MRSA from HCW to patient and vice versa was equal. In vitro studies have shown 
that the growth rate of certain CA-MRSA strains is faster compared to HA-MRSA strains 
suggesting that CA-MRSA may have an advantage towards colonization and therefore 
greater transmissibility . Model simulations using greater transmission parameters 

for CA-MRSA showed similar conclusions to the baseline model except for more rapid 
dynamics of CA-MRSA spread and emergence of co-colonization. Lastly, a deterministic 
model was used which compartmentalizes patients into homogeneous groups and therefore 
individual-level behavior was not addressed. Although stochastic individual-based models 
can simulate the heterogeneity of patients and HCW interactions, the increase in behavioral 
details provided by these models may result in greater difficulty in interpretation of findings. 

Colonization with different strains and species of bacteria is common [9]. This 
study quantified the emergence and spread of co-colonization with both CA-MRSA and 
HA-MRSA. The study demonstrated that the expanding community reservoir of CA-MRSA 
resulting in an increase infiux of CA-MRSA patients into the hospital setting coupled 
with prolonged LOS associated with severe CA-MRSA infections will rapidly lead to a 
predominance of patients who are colonized with both strains. The impact of these findings 
on patient outcomes and the potential for transfer of genetic information between these 
strains will require ongoing evaluation. 
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Appendix 

A reduced version of DAgata et al [t^ was developed, wliicli exhibits the same 
qualitative properties as the original, but which lacks the infective states. The patient 
population is split into three compartments: S - susceptible patients, C - patients colonized 
with CA-MRSA and H - patients colonized with HA-MRSA. We added the assumption 
that the total number of patients was conserved at size N = 400, allowing us to reduce 
by one dimension to the compartments C and H, by letting S = N — C — H — B. The 
model is analyzed without this assumption in an accompanying paper, and the results are 
qualitatively the same. 

We then extended the model, allowing patients to be concurrently co-colonized with 
CA-MRSA and HA-MRSA, which adds a third (or fourth if there is no conservation) 
compartment, B - patients colonized with "B"oth CA-MRSA and HA-MRSA. The rate of 
change of the size of the compartments due to the transmission of MRSA in the hospital is 
then described by the following system of nonlinear ordinary differential equations: 



^ = i6sS + 5cC + 6hH)Xc + iL^^siC + B) - ^1-^^^C{H + B) - {5c + ac)C 
at N i\ 

(1) 

^ = {6sS + 6cC + 5hH)\h + ^^^^^S{H + B) - ^1^Z^^£^H{C + B) - {5h + aH)H 
at N N 

(2) 

^ = i6sS + 6cC + 5nH)Xs + ^^^^^C{H + B) + ^1^^^H{C + B) - (6^ + as)B, 
at N N 

(3) 

with S = N — C — H — B. Parameter explanations and values are given in Table [TJ 



Assuming that all patients enter the hospital susceptible (Ac = = A^ = 0), 
there exists a disease free equilibrium (DFE) of the system, where all patients are 
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susceptible, S = N and (C, H, B) = (0, 0, 0). By linearizing the system and evaluating the 
eigenvalues of the Jacobian matrix, we find that the DFE exists for all parameters and is 
locally asymptotically stable if max{RO^ , RO'^} < 1, where RO^ = is the basic 

reproduction number for HA-MRSA and RO'" = ^l^^^c basic reproduction number 

for CA-MRSA. This means that if max{RQ^ , RO'^} < 1, both strains of MRSA will be 
extinguished over time. 

In addition to the DFE, there are two other analytically known equilibria. Eh and Ec, 
which describe states where one disease is endemic while the other absent. When E^ is 
stable, HA-MRSA will be endemic in the hospital and CA-MRSA will be extinguished over 
time. Therefore, the size of the compartments S and H will be positive and the size of 
compartments C and B will be zero. Symmetrically, when Ec is stable, CA-MRSA will be 
endemic and HA-MRSA will be extinguished over time. 

A fourth equilibrium, in which the size of every compartment is positive, does not have 
a known analytic form but is consistently found in numerical simulations. 



The invasion reproduction number [2^ 10 is a threshold parameter which determines 
the stability of the equilibria Eh and Ec- When 10^^ > 1, Eh is unstable, CA-MRSA 
invades, or becomes endemic, and both strains become prevalent in the hospital over time. 
Conversely, if JO^^ < 1, Eh is stable and only HA-MRSA will be endemic in the hospital 
over time. Symmetrically, when 10^^ > 1, Ec is unstable and both strains become endemic 
over time. Conversely, if JO^'^ < 1, Ec is stable and only CA-MRSA will be endemic in the 
hospital. 



A symmetric form is found for 10^^ , since the model is symmetric in C and H. 



(4) 
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Assuming fie — — Pen — Phc and ac — och = ckb, as we do in this paper, 

7o^« ^ — r- + - 1 

1 Kr H Kr 

When RO" > 1 and i?0^ > 1, /0^« > 1. 

Symmetric results hold for 10^"^ . Therefore, neither endemic equilibrium is stable when 
RO^ > 1 and i?0^ > 1, and there is never competitive exclusion. If patients continually 
enter the hospital colonized with MRSA, then it can never be completely extinguished. In 
an accompanying paper, the model is more thoroughly investigated mathematically. 
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Table 1: Parameter estimates for the transmission dynamics of community-acquired and 
hospital-acquired methicillin-resistant Staphylococcus aureus colonization (CA-MRSA and 
HA-MRSA). 



Parameter 


Symbol 


Baseline Value 


Source 


Total number of patients 


N 


400 




Percent of admissions per day 








Colonized CA-MRSA 


100 Ac 


3 


13,14 


Colonized HA-MRSA 


100 Xh 


7 


BI, 13, 14 


Length of stay 








Susceptible 


l/Ss 


5 days 


BI 


Colonized CA-MRSA 


l/Sc 


5 days 


BI 


Colonized HA-MRSA 


1/Sh 


7 days 


5 


Co-colonized 


1/6b 


7 days 




Hand-hygiene compliance efficacy (as %) 


100 r] 


50% 




Transmission rate per susceptible patient to 








Colonized CA-MRSA per colonized CA-MRSA 




0.4 per day 


11,12 


Colonized HA-MRSA per colonized HA-MRSA 




0.4 per day 


11,12 


Transmission rate per patient colonized with CA-MRSA to 








Co-colonized per colonized CA-MRSA 




0.4 per day 


11,12 


Transmission rate per patient colonized with HA-MRSA to 








Co-colonized per colonized HA-MRSA 


Phc 


0.4 per day 


11,12 


Decolonization rate per colonized patient 








per day per length of stay (as %) 








CA-MRSA 


100 ac 


0% 


15,16 


HA-MRSA 


100 a„ 


0% 




Co- Colonized 


100 as 


0% 





BI: data obtained from the Beth Israel Deaconess Medical Center 
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Fig. 1. — A compartment model describing the transmission dynamics of CA-MRSA and 
HA-MRSA in a 400-bed hospital. The arrows and parameter values correspond to entry and 
exit from the 4 compartments (^'-susceptible patients, C-patients colonized with CA-MRSA, 
if-patients colonized with HA-MRSA, and S-patients co-colonized with both strains). The 
percentages of patients admitted colonized with CA-MRSA, colonized with HA-MRSA, or 
colonized with both strains are expressed as lOOAc, IOOXh, and IOOAb, respectively. Dis- 
charge and death rates from the compartments are expressed as follows: Ss, Sc, Sh, and 
6b for susceptible patients, patients colonized with CA-MRSA, patients colonized with HA- 
MRSA. and patients co-colonized with both strains, respectively (with mean length of stays 
defined as 1/^5, l/Sc, 1/5^, and 1/Sb)- The colonization rates of susceptible patients to the 
CA-MRSA compartment is {l-r])f3c and to the HA-MRSA compartment is {l-r])f3H- The 
co-colonization rate from C to the co-colonized compartment (B) is (1 — r])(3cH and from H 
to i? is (1 — r])l3HCi where 100?] signifies the percentage of hand-hygiene compliance (where 
77 = corresponds to 0% compliance and r] — 1 corresponds to 100% compliance). The rates 
of decolonization of patients with CA-MRSA, HA-MRSA, or both strains are given by ac, 
atf. and an. resoectivelv. 
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Fig. 2. — Time evolution, over two years, of the percentage of patients colonized with HA- 
MRSA (dashed), CA-MRSA (dotted) or both (solid). Length of stay for HA-MRSA and 
both is 7 days, and 5 days for CA-MRSA. Initially the patient population consists of 90% 
susceptible patients, 3% of patients colonized with CA-MRSA, 7% of patients colonized with 
HA-MRSA, and 0% co-colonized. 
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Fig. 3. — Top - percentage of patients colonized with HA-MRSA (dashed), CA-MRSA 
(dotted), and both (sohd) after 2 years, as transmission {(3 = Pc = Ph = Pch = Phc) 
is increased. Bottom - i?0^ (dashed) and RQp (dotted) as transmission is increased. The 
invasion reproduction number 70^'* is plotted with the dash-dotted hne (see Appendix for 
description) . 
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Fig. 4 — Left - percentage of patients colonized with HA-MRSA (dashed), CA-MRSA 
(dotted), and both (sohd) after two years, as the percentage of patients entering the hospital 
already colonized with CA-MRSA (Ac) is increased. Right - the total percentage of patients 
colonized with MRSA after 2 years as Ac is increased. In both figures, IOOXh — 7% and 
lOOAs = 0%. 
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Fig. 5. — Left - percentage of patients colonized with HA-MRSA (dashed), CA-MRSA 
(dotted), and both (sohd) after 2 years, as the length of stay of patients colonized with 
CA-MRSA is increased. Right - the total percentage of patients colonized with MRSA after 
2 years as LOS is increased. LOS of patients with HA-MRSA is 7 days. LOS of patients 
colonized with both strains equals the greater of the LOS of patients colonized with CA- 
MRSA or HA-MRSA. 
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Fig. 6 — Top - percentage of patients colonized with HA-MRSA (dashed), CA-MRSA 
(dotted), and both (sohd) after 2 years, versus two interventions: left - decolonization (a) 
and right - hand- hygiene compliance {r]). Bottom - sum of percentages of patients colonized 
with either or both strains of MRSA. lOOAc = 3%, lOOAj^ = 7%, IOOAb = 0%. 



